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Beam position monitors (BPMs) have been widely used in all kinds of measurement systems, feedback sys- 
tems and other areas in particle accelerator field these days. The malfunction of a single BPM can cause serious 
consequences such as the failure of the orbit feedback and the transverse feedback. A troubleshooting has been 
made to prevent the defective BPMs from affecting the accuracy and stability of the storage ring in Shanghai 
Synchrotron Radiation Facility (SSRF). Different types of malfunctions have been successfully identified by 


using the idea of principal component analysis (PCA). 
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I. INTRODUCTION 
A. Various of Defective BPMs in SSRF 


A beam position monitor (BPM) system, is an essential di- 
agnostics in storage ring of a light source. The storage ring in 
SSRF is equipped with 140 BPMs located at 20 cells of the 
storage ring to monitor the beam dynamics [1]. The BPMs 
at the beam lines after the insertion devices (ID) or the bend- 
ing magnets are of great importance, because they also serve 
as the orbit feedback system to ensure stability of the elec- 
tron beams. BPMs of low resolution must be excluded from 
the feedback system for high stability of the beam. Correla- 
tion analysis based on principal component analysis (PCA) is 
used for disqualifying BPMs having low confidence or being 
considered as faulty or noisy. The BPM confidence levels in- 
cluded in the feedback system can also be used to estimate 
stability of the beam dynamics. Some BPMs can be used to 
do measurements other than the beam position, such as the 
(relative) beam current or life time. Therefore, an abnormal 
BPM should be found and treated. 

A typical BPM system consists of the probe (button-type 
or stripline-type), electronics (Libra Electronics/ Brilliance in 
SSRF) and transferring component (cables and such). Ever 
since the SSRF commissioning in 2009, the following BPM 
noise-causes have been found: 1) permanently damage to in- 
dividual probe or corresponding cable, 2) misaligned (posi- 
tion/angle) probes, 3) high-frequency vibrations, 4) electron- 
ics noise, and 5) others. 

A damaged probe or cable means totally useless of the sig- 
nals from the BPM, which should be ignored until its replace- 
ment or repair. A misalignment causes gain shift of the probe, 
hence the need of recalibration. The vibration and electronic 
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noise of a BPM result in increase of its measurement varia- 
tion and uncertainty. PCA has been an efficient and sufficient 
tool to locate the BPM malfunctions at SSRF. 


B. PCA 


PCA has been a useful statistical technique for finding in- 
terested patterns within mass data in various fields. It was in- 
troduced to the society of particle accelerators as a major con- 
stituent part of the model-independent analysis (MIA) theory 
by Irwin J and Wang C X in Dr. Wang’s doctoral work [2]. 
A series of work inspired by MIA have been done since then 
such as the optical parameters measurement [3-7] and per- 
formance estimate [8, 9]. The correlated BPM data matrix is 
regarded as a linear combination of different physical modes 
at different locations. A statistical analysis is then applied to 
study the beam dynamics without knowing the lattice model 
of accelerator. 

Assuming the following form of BPM output signals: 


b;(t) = alm, (t), (1) 


where, b;(t) is signal of the i BPM, m,(t) is the j" physical 
mode and a? is coefficient of the j™ mode. The Einstein sum- 
mation convention is adopted here to simplify the expression. 


It is found that the following physical modes are shared by 
all the BPMs in SSRF: 


1) A major contribution to the beam dynamics is the be- 
tatron oscillation. It has two degrees of freedom, so 
two physical modes should be sufficient to determine 
it: m1(t) = sin(wt + o) and ma(t) = cos(wt + ¢o), 
where ġo is a trivial initial phase that is shared by all 
BPMs. The corresponding coefficients a} and a? in 
Eq. (1) are the @ function and should be identical or at 
least relatively close to each other. 


2) The energy oscillation is visible. 
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3) Global signals are occasionally observed from the RFs, 
the power supplies and other parts of SSRF. 


Generally, variables of the physical phenomena caused by 
different physical mechanisms are considered independent, 
and independent variables are always linearly uncorrelated. 
Then, m,’s have formed an orthogonal basis which can fully 
describe the beam dynamics in its vector space. PCA is tool 
to find the basis. By using the singular value decomposi- 
tion (SVD), the BPM matrix B can be decomposed into three 
matrices [2]: 


Bmxn = ana og (2) 


where U and V are unitary square matrices and S is a diag- 
onal matrix. In MIA, U is called the temporal matrix, V the 
spatial matrix and S' the singular value matrix for the follow- 
ing reasons: 


e U is the matrix of the eigenvectors of the covariance 
matrix BB? and each column of U corresponds to the 
time series waveform of a specified mode; 


e V is the matrix of the eigenvectors of the covariance 
matrix BİB and each column of VÝ corresponds to 
the spatial distribution of a specified mode; 


e Each element of S is nonnegative and real, and it cor- 
responds to the variance of a specified mode. 


The decomposition can be expanded in a more visual way: 
i k ink 
b; = SKURU» (3) 


where, b, s, u and v are elements of the matrices B, S, U and 
VÝ. Eq. (3) implies that as independence is a stronger form 
of linear independence, the BPM data consist of a complete 
set of basis which corresponds to the independent physical 
modes suggested by PCA. Thus, Eq. (3) can be further ex- 
panded as: 


i _ (81,4 , 68,1 B2 i | 8,2 n, i n 
bj = 581UB 1j + Spaug, 205” + SyUntj (4) 
EN, i „EN , N i, N 
+ SENUENU; + SNUNY; , 


where, 8 denotes the two modes of the betatron oscillation, 
7 denotes the energy oscillation, EN denotes the electronics 
noise and N denotes all other noise. By examining the eigen- 
vectors and the singular values, different modes can be further 
analyzed. 


II. EXTRACTED MODES FROM SSRF DATA 
A. Betatron Oscillation 


The betatron oscillation is usually suppressed by the trans- 
verse feedback system and may not be quite significant. To 
evaluate the relation between the BPM signals and calculation 
model, a kicker was used in an experiment on November 1, 
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Fig. 1. (Color online) Waveform (a) and spectrum (b) of the betatron 
oscillation. 


2010, to invoke the transverse oscillation. Two components 
of the almost same singular values were extracted and one of 
the components is shown in Fig. 1(a). Fig. 1(b) shows corre- 
sponding spectrum of this component, indicating that it is the 
betatron oscillation in the frequency of the horizontal tune. 
The other component had the same characteristics except for 
a 7/2 phase shift. The two components can form a complete 
basis in the betatron oscillation vector space, and the 8 func- 
tion can be derived by using the spatial vectors. 

By comparing the measured and calculated 8 functions, the 
BPMs that did not work properly can be found (Fig. 2(a)). 
Moreover, the relative difference of each BPM (Fig. 2(b)) can 
be used to make a list of confidence levels of the BPMs. The 
No. 68 BPM is a reference BPM. With a single button probe 
as four identical input channels, the No. 68 BPM does not 
contain the beam dynamics information, theoretically. As 
shown in Fig. 2, BPM No. 68 does not reflect the betatron 
oscillation behavior. 


B. Energy Oscillation 


Another beam dynamics mode extracted by PCA is the en- 
ergy oscillation. It is of lower frequency, and the phase of 
oscillation tends to be invariant. In other words, the mode has 
the same behavior for all the BPMs within one resonance pe- 
riod. Thus, the “location free” component is one-dimensional 
in the vector space. A typical waveform of the energy oscilla- 
tion on July 22, 2010, was shown in Fig. 3(a). The spectrum 
of basis in Fig. 3(b) confirms that the mode is limited in the 
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Fig. 2. (Color online) Comparison of the measured and modeled 
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Fig. 5. (Color online) Waveform (a) and spectrum (b) of the 9" 


Fig. 3. (Color online) Waveform (a) and spectrum (b) of the energy i 
mode. 


oscillation. 
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Fig. 6. (color online) Spatial distribution of the 9" mode. 
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Fig. 7. (Color online) Comparison of the measured and calculated 
horizontal -functions (a), and the standard deviations (b). 


low frequency area. 

Following the procedures in Sec. II, and comparing the cal- 
culated and measured dispersion function, we found that the 
malfunctions corresponded to the spatial vector of the com- 
ponent (Fig. 4). 


C. Noise 


After an SVD, the modes were separated. Let’s take a close 
look at, for example, the mode corresponding to the 9" sin- 
gular value (referred as the 9"" mode) of the data on April 27, 
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Fig. 8. (Color online) Comparison of the measured and calculated 
vertical G-function (a), and the standard deviations (b). 


2011. The waveform in Fig. 5(a) of the mode does not seem 
to show much inspiring, but the spectrum in Fig. 5(b) demon- 
strates characteristics of the mode: it is a combined signal of 
the horizontal betatron oscillation and some 29 kHz electron- 
ics noise, and the vertical betatron oscillation was coupled 
into the mode as well. The spatial distribution of the mode 
gives the weights in all the BPMs. Fig. 6 gives a rough list of 
the BPMs which were suffering from this kind of noise and 
they are not suitable for extreme precise measurements. One 
sees that although BPM No. 68 is absolutely not related to the 
beam and gives only false signals, with its use of a relatively 
stable electronics. 


HI. USING PCA IN BPM TROUBLESHOOTING 


The BPM malfunctions mentioned in Sec.I do not behave 
alike in the PCA. Permanent damage to a cable or probe 
makes a total mismatch of measured physical parameters to 
their theoretical data. A probe misalignment causes shifted 
results of parameter measurement, but this is still adjustable. 
The vibration or electronics noise, however, may introduce 
some new modes which are shared by some of the BPMs. A 
series of BPM turn-by-turn data matrices were recorded on 
April 27, 2011. The results are used here to illustrate the 
BPM troubleshooting. 

The betatron oscillation was invoked by using a kicker 
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Fig. 9. (Color online) Comparison of the measured and calculated 
dispersion function (a), and the standard deviations (b). 


magnet so that the first two major principal components af- 
ter the decomposition of the BPM matrix were the betatron 
oscillation. The third mode was the energy oscillation. The 
BPMs that were suffering from the electronics noise (the 9" 
mode) had already been shown previously in Figs.5 and 6. 

The £ functions and dispersion functions were averaged, 
and then compared with the model in the first phase as shown 
in Figs. 7(a), 8(a) and 9(a). The relative differences between 
the measured parameters and the model ones are shown in 
Figs.7(b), 8(b) and 9(b) to illustrate the fitnesses of the BPMs. 
It is obvious that BPM No. 68 did not share the main beam 
dynamics and the corresponding (-functions and dispersion 
function (the spatial vectors) were all zeros which were not 
approximate to the mode values. 

Histograms of the measured horizontal, vertical ø- 
functions and the energy oscillations had been made after a 
series of measurements (Fig. 10). The BPMs with larger vari- 
ances of these parameters indicate that there might be a mis- 
match between the RF and the ADC sampling clock or the 
DDC local oscillator of the very electronics. 

The BPMs shown in Fig. 6 were believed to have serious 
electronics noise and this BPM list is not suitable for accurate 
measurements. The troubled BPMs shown in Fig. 4(a) were 
not necessarily useless in general. The ones with surprisingly 
small variances in Fig. 4(b), e.g., BPM number 123, might 
due to poor configurations of the electrodes, alignments or 
cable connections which contributed some offset to the results 


Nucl. Sci. Tech. 25, 020102 (2014) 


0.2 


(a) 


> 

> ha 

_ o 
T T 


= 
> 
a 


Measured fx deviation / Bx 


BPM ID 


Measured fy deviation / By 


BPM ID 


Measured y deviation / n 


40 60 80 
BPM ID 


Fig. 10. (Color online) Measurement variances of the horizontal 8- 
function (a), vertical $-function (b) and the dispersion function (c). 


and were hopefully removable or adjustable. 


For those unusable BPMs, not only the extracted Twiss pa- 
rameters were not close to the designed ones, but also the 
measured results were very unstable, i.e., had very large vari- 
ances. This was due to rapid changes of the local lattice 
or serious problems of the detector system. The lattice was 
presumably stable during the normal operations so that these 
BPMs were marked as damaged for future tests. 


All these measurements and statistics had been used to con- 
struct a list of confidence levels. The ones with higher confi- 
dence levels would be in higher priority in serious measure- 
ments or feedbacks. The misaligned ones would be calibrated 
once more to get rid of the gain shifts. The damaged ones 
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would be rechecked and cables would be replaced if the other 
parts of the probe would still work. 


IV. CONCLUSION 


By extracting the physical signals from the BPM data ma- 
trix, PCA has been proved to be a useful tool to separate vari- 
ous malfunctions in the BPM systems in SSRF. The machine 
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